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Abstract 

We show that in dilute boson-fermion mixtures with fermions in two internal 
states, even when the bare fermion-fermion interaction is repulsive, the ex- 
change of density fluctuations of the Bose condensate may lead to an effective 
fermion-fermion attraction, and thus to a Cooper instability in the s-wave 
channel. We give an analytical method to derive the associated T c in the 
limit where the phonon branch of the Bogoliubov excitation spectrum of the 
bosons is important. We find a T c of the same order as for a pure Fermi gas 
with bare attraction. 
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The interest in the effective interaction between fermions in boson-fermion mixtures is 
not new. Already in the 1960 's Pines |I| suggested that the effective interaction between 
3 He atoms in solution in superffuid 4 He could be attractive due to the exchange of density 
fluctuations of bosonic background. This attraction was then observed experimentally by 
Edwards et al. [0] and Anderson et al. ||. From the experimental data Bardeen, Baym and 
Pines @] estimated the expected critical temperature for Cooper pairing both in the s-wave 
and p-wave channel. The situation is similar to that in ordinary superconductors where the 
effective attraction between electrons is well established to be caused by the exchange of 
lattice phonons ||. 

The renewed interest in the issue stems from the recent availability of trapped atomic 
gases at the ultralow temperatures required for quantum phenomena to be relevant. After 
Bose-Einstein condensation (BEC) was observed, huge steps have been made also in the 
cooling of fermions |J ||10|| . One way to obtain the cooling |10] has been to mix 6 Li 



atoms (fermions) with 7 Li atoms (bosons), and to proceed with the standard technique of 
evaporative cooling on the latter so as to cool the former 'simpathetically', i.e. indirectly, 
by simple thermal contact. At the end of the process a stable nearly pure BEC of 7 Li on 
top of a cloud of degenerate 6 Li has been observed. It is therefore timely to calculate the 
effective interaction between fermions in a dilute boson-fermion mixture when fluctuations 
of the BEC are taken into account, and to study the conquences on Cooper pairing. 

In two recent works Heiselberg, et al. |Tl| and Bijlsma, et al. |12j] found the same expres- 



sion for the boson-induced contribution to the interaction, but proceeded in different ways 
to analyze the effects on s-wave pairing. The equivalence of the two approaches in dealing 
with the latter problem is not obvious since the results were given in very different forms, 
one analytical [Tl[ and the other numerical [12]. Moreover, their applicability is, for ref. [II 



limited to the case of an already attractive bare fermion-fermion interaction, and for ref. (12 
limited by the lack of transparency and physical understanding of numerical calculations. In 
this paper a third way to treat the problem is developed. The idea is to give an approximate, 
but under usual circumstances accurate, analytical solution to the numerical approach of 



ref. [0. This applies also to the case of a bare fermion-fermion repulsion, and we prove 



that, as conjectured in []I2J and as one expects, even if the bare interaction is repulsive the 
effective one may be attractive when the boson-induced interaction is taken into account. 
In addition to that we find an elegant and simple expression for the associated BCS critical 
temperature, which contrary to the case of 3 He- 4 He mixtures is fully determined, also in 
the prefactor. Since our approach is valid for bare fermion-fermion attraction as well, we 
also show that it gives, within its range of applicability, the same results as were found in 



ref. |p-l|| , and thus prove the consistency of the three methods. 

When the gases of the mixture are dilute all interactions can be described by one pa- 
rameter each: the appropriate scattering length [|13"| , |14|| . In a mixture of fermions in two 
different internal states 1 and 2 (with the same mass), and bosons in one internal state only, 
the scattering lengths a 12 , cibb, a Bi and a,B2, characterize all the relevant interactions, since 
the Pauli principle allows us to ignore interactions between identical fermions at the tem- 
peratures and densities of interest here. In the following we shall suppose, without loss of 
generality, that asi = a B2 = a BF , and rename ai 2 = a FF . Often in the text pseudopotentials 
will be used instead of scattering lengths. They are defined as follows: U FF = 4:irh 2 a FF /m F , 
Ubb = 4:Kh 2 aBB/ m B, and Ubf = 4:^h 2 a BF /m BF , where m BF = 2m B m f / \m B + m F ) is 
twice the reduced mass for a boson with mass m# and a fermion with mass m F . Since we 
are interested in s-wave pairing we shall suppose that the densities of 1 and 2 fermions are 

the same, n\ = n 2 = n F , and therefore also their Fermi energies and momenta, e F and hk F . 
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Moreover we shall introduce in the theory small parameters of the type k F a and n B a, where 
by a we indicated a generic length of the order of a scattering length and by n# the boson 
density. 

In a mixture the effective interaction between a spin 1 and a spin 2 fermion is the sum of 
the direct one U FF , the one arising from polarization of the bosonic medium U F bf, and that 
due to polarization of the fermionic medium itself U FFF . We shall show that under suitable 
conditions the first two interactions can be of the same order, and it is therefore essential 
to consider both of them to predict whether the system undergoes Cooper pairing or not, 



and if so what is the critical temperature T c . The third one is instead of order a(kpa) and 
is important for the renomalization of the prefactor in the expression for the BCS critical 
temperature. The way to deal with Ufff has been shown in ref. ]15[ (see also There 



it was found that it causes a decrease in the prefactor of T c by a factor (4e) 1//3 , and we take 
this for granted here. 



The boson-induced interaction between fermions introduced in refs. [11,12 



is 



U FBF (q., w) = U BF x(q_, w), (1) 

where the boson density-density response function in the Bogoliubov approximation is given 
by 

X <l w =7F^ of o i o — FT — v 2 

(huj) 2 - e° q (e° q + 2n B U BB ) 

e q = h 2 q 2 /2m B and q — |q|. Since we are considering dilute Bose and Fermi gases, we have 
neglected the renormalization of x(q, due to the presence of the fermions, this is correct 
to the lowest order in the gas parameter. We used the zero temperature response function 
because the BCS critical temperature is much smaller than the boson condensation one if 
n B > n F . 

The interaction in Eq. ([I]) provides an attraction between two particles at the Fermi 
level, since in that case u = and Ufbf{q,Q) < 0. We remind the reader that in uniform 



systems U BB > is required for stability of the mixture [|T(| . 

At this point one has to analyze two possibilities. The bare fermion-fermion interaction 
can either be attractive Uff < or repulsive Uff > 0. In the former case the gas would 
undergo pairing even in the absence of bosons at the critical temperature T c . In the presence 
of bosons the direct and induced contributions add up constructively to a stronger attractive 
interaction, and the BCS critical temperature rises. This possibility was studied in detail in 



ref. ||11|| . If Uff > instead, the Fermi gas in not unstable to pairing without the bosons. 
When the bosons are added however, if the attractive boson-induced interaction at the Fermi 
surface is stronger than the bare repulsion, the total effective interaction Uff + Ufbf(q, 0) 



is attractive, and the gas becomes unstable to pairing. One can immediately recognize 
a mechanism at work completely analogous to that of electrons in superconductors. Just 
as in superconductors the induced interaction depends on the energy exchanged, and is 
attractive only in a band centered about the Fermi surface 0. Obtaining a solution for 
T c for arbitrary densities is complicated. As is well known however PJT^] if vp <C s, where 



s = {j^bU bb I ttib) ^ is the sound velocity in the Bosegas, retardation effects can be neglected 
and to can be set to zero always. The induced interaction is thus attractive in the whole 
Fermi sphere, and we seek a solution to the problem under this assumption. 
The boson induced interaction when vp <C s is then 

Ufbf ® = -fh ' i I (* To Vr ^ 

Ubb 1 + {nq/2m B sy 

Notice that if m# > mp, since the typical momentum exchanged in an interaction is fiq ~ 
mpVF <C 2rriBS, we expect Ufbf(q) — ~U bf /Ubb, i-e. a constant independent of q. 

We now want to properly take into account both bare and induced interactions to deter- 
mine T c . According to Emery |T7| if the fermions with opposite spins interact via a potential 
U(r), and if tan<5 (&F) > 0, where S^kp) is the associated / = phase shift evaluated at 
the Fermi wavenumber, the system undergoes s-wave pairing at the critical temperature 

k B T c = 1 (-) 7/ \p e -V2t*»*>(M ( 4 ) 



7r \ej 

In writing Eq. we have already included the correction to the prefactor due to the 
polarization of the fermions. 

In the case of a pure two-species Fermi gas with bare attraction, U(r) is the bare potential 
(with associated scattering length and pseudopotential app < and Uff respectively), and 
since by assumption k F \aFp\ 1, then tsna6 (kp) ~ —kpapF, and Eq. (H) reduces to the 



well know formula |T5j 
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k B T c = 1 (-) e F e n/2kpaFF . (5) 
In a mixture on the other hand U(r) changes. This approach for finding the effects of 



the boson- induced interaction has been used by Stoof and co-workers |12|. To apply Eq. 



to our case one has to take the Fourier transform into real space coordinates of (|3]) 

UFBF(r) = --^- T . 1 -e-^, (6) 
AitUbbZb r 

where = ^/2m^s is the boson coherence length. The 1/r divergence at r = is artificial, 
since the potential must in any case be cut off at a distance r of the order of a scattering 
length. Ufbf{t) is a Yukawa potential with range = asB (167T nsa^ B ) _1//2 which is much 
greater than orb if the gas is dilute. 

The fermions interact both via £/bare(j) and Ufbf{t) so that the total interaction po- 
tential is given by U tot {r) = t/barc( r ) + Ufbf(j')- The aim is then to calculate the s-wave 
phase shift due to the total potential U tot {r). In principle this is a difficult problem, since 
one should solve the radial Schrodinger equation 

/ d 2 nip 



, (Ir2 +-^U tot (r)-k 2 ju™(r;k) = 0, (7) 

where, however, the bare potential is not known, as only the scattering length app is the 
measured quantity. Fortunately it is not necessary to know also the details in r of £7bare- 
Eq. (0) can be largely simplified by noticing that £/b are and Ufbf act on two different length 
scales. The first one from r = to r ~ r$ and the second one from r ~ r$ to r ~ ^> tq. 
One can then just solve Eq. (0) for r > r , with U tot {r) = UpBF( r ), and introduce a 
boundary condition on MQ 0t (r; k) at r = r , which accounts for the phase shift due to the 
bare potential. This is due to the fact that by the time the wave function reaches the region 
of distances where UpBF{ r ) is relevant, the bare potential has stopped acting and the wave 
function has recovered its sinusoidal form with phase shift <5g are (/c). 

Using these replacements Eq. ([7]) has been solved numerically by Bijlsma et al. (12| for 
various sample elements and densities, and the results for #o ot (^) are in their publication. 
Whenever tan<5Q 0t '{kp) > the critical temperature is found by replacing the value obtained 
into Eq. (|). 

We shall now show that if kp^s C 1 (a condition automatically satisfied if mg > rap, 
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since by assumption v F <C s) tan <5o ot (/cp) can in fact be found analytically, and the final 
result is remarkably simple. 

It is well known, see for instance ||18|| , that given any two potentials U^(r) and U^(r) 



for the interaction of two particles with reduced mass m re d, and given the solutions with 
wavenumber k to the corresponding radial Schrodinger equations: Ui(r;k) and u\ 2 \r;k) } 
the associated phase shifts S^^k) and $\ 2 \k) are related by 

tan<y, (1) (Jfe) -tan& (2) (Jfe) = -k^^ H uf\r-k) 

n Jo 

x [[/« (r) - U {2) (r)] u\ 1} (r; k) dr. (8) 



This formula can be applied to our case by letting U^'(r) = U tot (r) and U^ 2 \r) = U^^r), 
m re d = rnp/2, and k = k F . For the I = channel then 

tan<5* ot (A; F ) = tanS^ re {k F ) - kp—j- (9) 



oo 

x / u h ave (r;k F )U FBF (r)ul ot (r;k F ), 



o 

and since k F \a FF \ <C 1, the bare phase shift is 

t&n5o are (k F ) ~ -k F a FF . 

We can now use the special form of our potentials. Because U F b F (t ) is zero for r < r , 
the relevant lower limit of the integral in Eq. (|9|) is really r . But for r > r we can replace 
tig are by its asymptotic value for r — > oo since the bare potential has by then decayed, thus 
MQ are (r; k F ) ~ kp l [sm{k F r) — k F a FF cos(k F r)}. If in addition k F ^ B <C 1, we can approximate 
Mg are (r; k F ) ~ r — a^p. 

Later we shall prove that if the gases are dilute and for typical values of the parameters, 
the potential U F B F {r) is shallow, i.e. not strong enough to form 'bound states'. Tjinhus we 
can to a first approximation (Born) also let «Q 0t (r; k) = iiQ are (r; k). This yields 

tan<$o ot (£; F ) = -k F a FF - A, (10) 

with 



m F k F U BF 



/•CO 1 

/ dr(r -a FF f-e- r/iB (11) 
J rn r 



{e- r "°(l + f ) -2a 



k F a FF U BF , _ 



UbbU ff 

~2 f°° , 1 -xl 

+ a FF / ax— e f 



= ?V£b, and a FF = a FF /£ B - 

The boson coherence length is in general much larger than the cut-off r which in turn 
is of order | a FF \ . Therefore f and a FF can be set equal to zero in the first two terms in the 
curly brackets. The last integral is dominated by the logarithmic divergence and it goes like 
~ a FF Info, but since a FF ~ f o <C 1 also the last term can be set to zero. The final result 
is therefore simply 

tanC(M = ~k F a FF (l - -% ) . (12) 

\ UbbUff J 

In order for the system to condense we need to have tan SQ 0t (k F ) > 0. This is always the 



case if a FF < 0, and we recover the result of ref. JTTJ] , in the limit k F ^ B 1. Notice that 



when k F ^B ~ 1 we cannot expand u^ re as we did anymore and we expect corrections to our 
result. However, the boson-induced attraction is maximized in the limit kp^B "C 1 since it 
is the phonon branch of the Bogoliubov spectrum that provides most of the attraction, and 
that is why it is satisfactory for the time being to consider only this limit. 

Moreover, what we found extends what is reported in ref. |TI[], since we can also see that 
when the bare interaction is repulsive, if U bf /UbbUff > 1 a system which would normally 
be stable undergoes a BCS transition due to the boson induced interaction. 



Our results are also consistent with those in ref. |12j since they were found using a 
limiting procedure of the same approach, but the values of the phase shifts are here given 
analytically, allowing a deeper physical understanding of their significance, and showing a 
simplicity which was hidden by the numerics. To be sure we also have solved the problem 
numerically and found the solutions given in ref. [TIJ. We have thus checked that they 
coincide with the present results in the limit kp^B *C 1. 

In all cases in which A is negative the critical temperature is given by 
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kBTc= l\e) eFC ' (13) 
where 

^'M 1 "^)' (14) 

and JV(O) = m F k F / (2vr 2 /i 2 ). This shows, as we already anticipated, that the critical tem- 
perature in the case of boson-induced pairing is of the same order as that of pairing with 
bare fermion-fermion attraction, and can be very large if U BF /UbbUff ^> 1- The condition 
can be achieved by a suitable choice of elements. At this time it is difficult to suggest an 
appropriate choice since the scattering lengths for collisions between bosonic and fermionic 
atoms are mostly being studied at the time of writing, and precisely in view of the present 
developments to which this work is a contribution. We also point out that an interesting 
consequence of ( |T3| ) and (|nj) is that, so long as k F ^s <C 1, the new critical temperature is 
independent of the boson density. The strong density dependence in the plots of ref. [0 is 
explained by the fact that the regime k F ^B ~ is also probed. 

In the derivation above we have used the Born approximation. This is valid only if the 
potential is sufficiently shallow. But, as we see from Eq. (§), the depth of the induced 
potential depends linearly on through the square of the coherence length at the denomi- 
nator. At low boson densities the potential is very weak and all the considerations above are 
certainly valid. But when tib is large enough for the induced potential to be able to host a 
'bound state' they fail, and the new scattering length and critical temperature may depend 
dramatically 011% as the calculations in ref. |I2| indicate. This regime is interesting and 
worth studing in greater detail, but is beyond the scope of this work. 

By a simple argument we may estimate the highest boson density allowed for our model 
to apply. The typical potential energy of a particle confined in the potential (|6|) is, apart 
from the sign, E p ~ U bf /^kUbb^,%i and the kinetic one E k ~ h 2 /m F ^ B . To be safe then we 
need to require E p /Ek <C 1, which implies 

n » aBB<< (16^{u^) ■ (15) 
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Thus if the scattering lengths and the masses are approximately the same, as in typical 
conditions, it is enough to require that the boson gas is dilute for the Born approximation to 
apply, but specific checks may be necessary especially for large Ubf- Recall that the boson 
density cannot be too low though if the condition Icf^e <C 1 also has to be fulfilled. 

In conclusion we have shown that the boson-induced interaction in a boson-fermion 
mixture can cause s-wave Cooper pairing in a Fermi gas with bare repulsion. We have 
calculated the associated critical temperature in the limit <C 1 where the highly efficient 
attraction due to the phonon branch of the Bogoliubov spectrum is important and found 
that T c is of the same order as for s-wave pairing in a gas with bare attraction. 

Finally, I express my gratitude to C. J. Pethick from Nordita and H. Smith from the 
0rsted laboratory in Copenhagen, were a large part of the work was carried out. I would 
also like to thank P. F. Bortignon and R. A. Broglia for support at the time of writing. 
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